Symmetry Protected Topological Order and Spin Susceptibility in Superfluid 3 He-B 
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We here demonstrate that the superfluid 3 He-B under a magnetic field in a particular direction 
stays topological due to a discrete symmetry, that is, in a symmetry protected topological order. 
Due to the symmetry protected topological order, helical surface Majorana fermions in the B phase 
remain gapless and their Ising spin character persists. We unveil that the competition between the 
Zeeman magnetic field and dipole interaction involves anomalous quantum phase transition where 
topological phase transition takes place together with spontaneous breaking of symmetry. Based 
on the quasiclassical theory, we illustrate that the phase transition is accompanied by anisotropic 
quantum criticality of spin susceptibilities on the surface, which is detectable in NMR experiments. 
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Introduction. — Superfluid 3 He-B is one of the most 
concrete examples of time-reversal invariant topological 
superfluids (TSFs) 0,0, where the ground state wave 
function supports a nontrivial bulk topological invariant 
in three spatial dimensions [H-IH . As a consequence of the 
bulk-edge correspondence, helical Majorana fermions live 
on its specular surface, and their self conjugate property 
gives rise to an Ising-likc anisotropy of spin susceptibil- 
ity Recently the surface Majorana cone has been 
detected in experiments fl7j ]. 

Since the topological superfluidity in 3 Hc-B, which is 
categorized to class Dili [3j|, is ensured by time-reversal 
invariance, it is sometimes stated that any time-reversal 
breaking such as a finite magnetic field immediately 
wipes out the topological nature. Indeed, in the presence 
of a strong magnetic field, the Majorana Ising spin ceases 
to exist [181 ]. However, as argued in the Letter, a more 
careful consideration on the basis of microscopic calcula- 
tions and symmetry of the system points to a different 
conclusion. It is worth mentioning that the robustness 
of a topological phase transition for a topological insula- 
tor against the time-reversal breaking is proposed in an 
extended Kane-Mele model, [lit . 

In this Letter, we show that 3 He-B under a magnetic 
field in a particular direction stays topological as a sym- 
metry protected topological order. In spite of the time- 
reversal breaking due to the magnetic field, the topo- 
logical property is retained by a hidden Zi symmetry 
that is obtained by a combination of time-reversal and an 
SO(3)x,+s rotation. Because this Z^ symmetry restores 
a chiral symmetry of the microscopic Hamiltonian, heli- 
cal surface Majorana fermions in the B-phase remain to 
be gapless, and the Ising spin character of the Majorana 
fermions persists unless the discrete symmetry is spon- 
taneously broken. Finally, we come to the conclusion 
that at a critical Zeeman field H* the system undergoes 
anomalous quantum phase transition in the sense that 
topological phase transition takes place together with 
spontaneous breaking of the Z^ symmetry (Fig. [Ha)). 



The phase transition is described by a pair of an or- 
der parameter l z and a topological number w, which are 
defined in Fig. QJb) and in the text below Eq. ([4]) . Con- 
ventionally, topological phase transitions are accompa- 
nied by creation or destruction of gapless surface states 
and gap closing of bulk quasi-particle spectra. While 
the anomalous quantum phase transition in the above 
involves the destruction of gapless surface states, it does 
not show the bulk gap closing. Instead, there appears a 
long range order due to symmetry breaking. Moreover, it 
takes place between two conceptional different quantum 
orders: At a critical magnetic field, the system undergoes 
a transition from topologically ordered state (w ^ 0) to 
conventionally ordered one (£ z ^ 0). 

The symmetry protected topological order is closely 
associated with the order parameter manifold of 3 He- 
B. Ignoring the dipole interaction and a Zeeman field, 
the bulk 3 He-B spontaneously reduces the symmetry 
SO(3)z,xSO(3) s xU(l) to SO(3)i, + s 0. The gap func- 
tion is A(k,r)~icr ■ d(k,r)a y , where d^ = d^ v k v and, 



(1) 



The broken symmetry SO(3)i,_s, the relative rota- 
tion between spin and orbital spaces, is described by 
R m (ii,if) with the rotation axis n and the angle ip. 
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FIG. 1: (Color online) (a) Schematic phase diagram of 3 He-B 
under a parallel magnetic field, where H* involves the topo- 
logical phase transition with spontaneous symmetry breaking, 
(b) Relation between t z and the the orientation of d(0, 0, k z ) 
for an arbitrary (n, ip) at the surface. 
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Here the repeated Greek indices imply the sum (/it, v, rj = 
x,y,z) and denotes the Pauli matrices in the spin 
space. In general, the dipole interaction acting as a small 
perturbation chooses a particular state of (h, ip). 

To quantitatively determine t z {n, ip) and H* , we here 
utilize the quasiclassical theory which takes account of 
dipole interaction and Zeeman energy on equal footing. 
In a slab geometry, the finite H* results from the compe- 
tition between dipole and magnetic energies, where the 
former (latter) favors the symmetry protected topologi- 
cal (a non-topological) order. It is found that since the 
topological order protects the Majorana Ising spins, the 
phase transition is accompanied by anomalous critical 
behaviors of spin susceptibilities on the surface. 

Surface bound states. — Let us start with the mean- 
field Hamiltonian density in the Nambu representation, 



K(n,r 2 ) 



e(n,r 2 ) A(n,r 2 ) 
-A*(n,r 2 ) -e*(n,r 2 ) 



Vz6(r 12 ).(2) 



In this paper, we set H = 
sists of e(r 1 ,r 2 ) = S(r 12 )(- 



&B = 1. Equation ^ con- 
-V 2 /2Af - E F ) and the Zee- 



man energy Vz = — fi n H IJt di&g(a l j,, — c^), where M, Ep = 
k F /2M, and fj, n are the mass, Fermi energy, and mag- 
netic moment of 3 He atoms. The pair potential for 
3 He-B, A(fc,r) = / dr 12 e~ tkri2 A(n, r 2 ) with Eq. (TJ, 
is simplified to A(fc, r) = U(n, p)Ao(k, r)U T (h, ip) with 
U(n,p) eS\J(2) and A (k,r) = ia IM a y A IM (r)k fJ ,. 

Wc first diagonalizc Eq. ((2]) as / dr 2 H{ri,r 2 )<pE(r 2 ) = 
E<Pe{ti), that is, the Bogoliubov-de Gennes (BdG) 
equation, where E and ip e describe the energy and wave- 
function of quasiparticles. This is solved within the An- 
dreev approximation V 2 — >ivp • V and the uniform pair 
potential A M (r) = Ao, where vp = kvp is the Fermi ve- 
locity. In this work, we consider the B-phasc sandwiched 
by two specular walls which are normal to the z-axis. 
For ff = 0, the dispersion of the surface Andreev bound 
state (SABS) is given by E (k») = ±^-\k»\ with fen be- 



ing the momentum in the xy-plane j8lll0l The cor- 
responding wave functions are expressed as (r) oc 
e ife r r|| e -z/$ gin ( J k 2 _ k^z)U(n, <p)^\ where (±) cor- 
respond to the positive and negative energy states and wc 



set U = 



Ediag({7, {/*). Wc also introduce <&j^ =2^*1 ^* = 

(l,-ie i<t>k ) -e i ' t>k ,-i)' T with (f> k = tan.- 1 (k y /k x ) and 
being the Pauli matrices in the Nambu space. 

For a finite H , the dispersion of the SABS is obtained 
from the linear combination, ip^ = a+(£>Q ^ +a_y>g fe ' , as 



E{k\\) = ± v /[So(fc||)] 2 + [Pnff4(n,^)] : 



(3) 



where i v (h,(p)=h fl R flu (n,cp) with h^ = H^/H. 

Symmetry protected topological phase — As we showed 
in Eq. ([3]), if i z = 0, the SABS remains gaplcss even in 
the presence of a magnetic field. From topological point 



of view, however, this seems to be a puzzle: Because 
the magnetic field breaks the time-reversal invariance, 
topological protection as a time-reversal invariant TSF 
does not work any more. Nevertheless, no gap opens in 
the SABS if the magnetic field satisfies £ z = 0. 

First, one should notice that £ itself could be affected 
by a Zeeman magnetic field. Therefore, an immediate 
solution for this puzzle might be that if one applies a 
magnetic field, i changes so as l z ^ 0. However, as is 
shown below, if the Zeeman field is parallel to the xy- 
plane (say, along x), this is not the case: There exists a 
symmetry that ensures £ z =0. Interestingly, we find that 
this symmetry resolves the puzzle above at the same time, 
by providing another topological protection of the SABS. 

Let us first consider symmetry of the system. Among 
SO(3)x, + s rotations under which microscopic interac- 
tions of 3 He atoms are invariant, the slab geometry con- 
sidered here preserves its subgroup SO(2)x, + s rotation 
U (0) in the xy-plane. The point is that while the Zeeman 
term along x explicitly breaks both of the time-reversal 
symmetry and the SO(2)x+s rotation symmetry above, 
it does not break a combination of them. In fact, in this 
case, the flipped magnetic field by time-reversal T is re- 
covered by the 7r-rotation in the xy-plane. Therefore, the 
microscopic Hamiltonian of 3 He atoms is invariant under 
the discrete symmetry given by 1~U(tt). Wc notice here 
that £ z mentioned above is transformed nontrivially as 
l z — > — £ z under this symmetry. Therefore £ z is an order 
parameter of the discrete symmetry, and it should be zero 
unless the discrete symmetry is spontaneously broken. 

Remarkably, one can introduce a topological invariant 
if the discrete symmetry is not spontaneously broken. 
In that case, the BdG Hamiltonian ([2]) in the momen- 
tum space is manifestly invariant under the discrete sym- 
metry, n{k xi k v , ~k z ) = Tu(ir)n(k x ,k y , k^u- 1 ^- 1 , 

where T = ia y K is time-reversal with complex conju- 
gate operator K and U (it) = icr z T z is the 7r-rotation. 
Therefore, combining it with the particle-hole symme- 
try of the BdG Hamiltonian, CH(k)tf = -H*(-k) with 
C = t x K, one obtains the relation VH (k T , k y , fc^)r -1 = 
~ k Xl — k y , k z ) with T_ = cr x T y . On the k z axis, this 
reduces to the so-called chiral symmetry 



{r,7£(0,0,fc 2 )} = 0. 



(4) 



Thus, following Refs. [2f| |27[, one can introduce the 
following one-dimensional (ID) winding number w = 
— r*L dk z ti\ T'H ~ 1 dk^'H] I . __, which are evaluated as 

'-± 7T L ■/ OO ~ I /C 1 1 — U 

w = 2 for fi n H < Ep (A± > 0). Therefore, the system 
is topologically non-trivial, and the bulk-edge correspon- 
dence implies that the SABS satisfies -E(fcj|)=0 at k\\ =0 
even in the presence of the magnetic field. It should be 
noted here that one needs the discrete symmetry spe- 
cific to this system in order to define w. Therefore a 
topological phase realized here is a symmetry protected 
topological order [2(J 21 1. 
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Majorana Ising Spin. — In the absence of the mag- 
netic field, it have been known that helical Majorana 
fermions in 3 He-B have Ising-like spin density Now 
we show that the symmetry protected topological order 
discussed above retains the Ising spin character of Majo- 
rana fermions. 

To see this, we use a general symmetry property of the 
SABS. Let us consider the low energy limit where only 
the zero energy SABSs at feii =0 contribute. According 
to the index theorem of Ref. [27| , Eq. (J4j) infers that the 
zero energy SABSs are eigenstates of T, and the relation 
w = n_ — n + holds on the surface of superfluid 3 He-B, 
where n± is the number of the zero energy SABSs with 
the eigenvalue T = ±1. In the present case, w = ri- =2 
and 7i + =0, and thus the SABS V^Lo satisfies ^ i P < k^=o = 
— VfeLo ( a = lj2). Using the particle-hole symmetry, one 



,( a )* 



(a) 



at the same 



can also put the relation T x cp k _q 
time. From these two relations, yj^Lo nas a generic form 

as VfejLo = K (q) - ^ {a) *,& } *, -i£ (a) ] T with a function ^ . 
Ignoring non-zero energy modes, the quantized field \f = 
[^Ai,^\,^\] T is expanded as #(z) = E a =i,2 ^k"=o7 (o) 
with real and from the general form of tp k -q, onc 

obtains iipf = — f/jj, which is a general consequence of our 
symmetry protected topological order. 

Now following Ref. 0, |T(| , one can show that the last 

/,t 



where denotes the Fermi liquid corrections obtained 
as i^(fc, r) = J2i Af(P e (k, k')g^(k, r; iWn)) k < tUn - P e is tnc 
Legendrc polynomials and (■ • ■ ) k ^ = TX)| w „| <Ec / £ 
with a cutoff E c . The coefficient A| = F//(1 + Ffj{U + 
1)) is parameterized with the antisymmetric Fermi liquid 
parameters, F a = -0.695 and Ff = -0.5 0. 

The gap equation is obtained as A ab {k, r) = 

(V r o c 6 d (fe,fe / )[* 'M C7 '»/A»(* / > r ;* a; n)]cd)fe' )U , nJ with a,&,c,d=t 
,4- At the low pressure limit, the pair interaction of 
3 He atoms is described as V£ b (k,k ) =3\A\k fl k'6 ac 6 b d — 
Q^u{k 1 k'){a^) ac {a v )bd- The first term arises from 
the p-wave interaction with SO(3)s x SO(3)x, x U(l) 
and the second term is the dipole interaction, where 
Q^{k,k') is obtained from Q M „(fc, fc') = .9D-R / r~ 3 (<5 M „- 
3f t tf v )e~ l ( k ~ k }' r dr with kfzkkp. The factor R includes 
the contributions of high energy quasiparticles [2^ | . The 
dipole interaction can be expressed in terms of the par- 
tial wave series (jp-, /-, and higher waves). However, 
since the pairing interaction between 3 He atoms is dom- 
inated by the SO(3) s X SO(3)x xU(l) channel and the 
dipole interaction can be regarded as a small pertur- 
bation, we take account of only the p-wave contribu- 
tion of Qn V (k,k'). Then, the gap equation for A M (r) is 



given by A„(r) = 
3A D [<5,^(fc, ) / J) )^_ 



relation, iij)^ = —ip], yields the Ising character of the limit with A^^A 

reproduces cos ip = 



"I 

SABSs: It is shown that among the local density op- 
erator and the spin density operators, which are given by 

P=\[^i^a-^a^l\ &ndS^ = \[^l(a^) ab ^ b -^ a {aJ) ab ^l\, 

respectively, only S z is nonzero while the other compo- 
nents are identically zero. So, in the low energy limit, the 
SABSs do not contribute to the local density fluctuation, 
and its local spin density is Ising-likc. Here note that 
we only use a general property of the chiral symmetry, 
thus the Ising character is a direct consequence of our 
symmetry protected topological order. 

Quasiclassical Eilenberger theory. — Let us now micro- 
scopically determine l z and H* . For this purpose, we 
here utilize the quasiclassical Eilenberger theory, which 
provides a quantitative theory for superfluid 3 He at low 
pressures [28[ . This is based on the quasiclassical Green's 
functions g = g(k,r;iuj n ) with the Matsubara frequency 
u> n = (2n + 1)ttT (n€Z) and the 2x2 unit matrix ctq 



S^R^(h lV >){(3\A\^A D )(k u f n ) Lun - 

<( fc X f)r)k J}> With the di - 

Vt the th 

Aj_, the gap equation 



mensionless factor An = ff<?D-R- At the thermodynamic 



=An 

1A 
4 A 



and A 
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The Eilenberger and gap equations with Eq. (fTJ) pro- 
vide self-consistent equations for g under a fixed (n, ip). 
The Eilenberger equation with g 2 = —tt 2 t is numeri- 
cally solved with the Riccati parameterization in the sys- 
tem that two specular walls normal to z are situated at 
2 = and z = 20£. The numerical procedure is same 
as that in Refs. [3, We solve the gap equation 



with IAI-^ttTcoE 



|w„|<b c 



\UJ„ 



and A D /A 2 =2xlO" 
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o-o9o + 
ia y fl+ia y o-^fl 



(5) 



(2xl0 -5 ), where we set the cutoff E c = 2QttT c0 . Note that 
since the ratio Ad / A 2 is associated with the distortion of 
A M in the thermodynamic limit (32j . it is independent of 
the energy cutoff. The coherence length £ = v-p/ irT c0 is 
estimated about 80nm at the zero pressure of 3 Hc-B and 
T c o ~ lmK is the critical temperature of the bulk B-phase 
under H — 0. In a slab geometry, n may be assumed to 
be spatially uniform, because the length scale of the spa- 
tial variation of n is macroscopically large [22[ , compared 
with the thickness of the sample 20£« 1.6/jm in typical 



The evolution is governed by the Eilenberger equation 
[ico n T z — S_(k,r),g] = —iv-p ■ Vg. The 4x4 matrix S_ 
consists of the Zeeman energy Vz and the self-energies, 



experiments [34|, [35J. Then, the stable configuration of 



S(fc,r) 



f 



;r z V z 



v (1 (k,r)o- fl A(fc, r) 
At(-fe,r) p;(k,r)a; 



(6) 



(to, ip) is determined by minimizing the thermodynamic 
potential 5Q[g] whose explicit form is same as that given 
by Vorontsov and Sauls in Ref. j2f|. Note that for the 
thickness 20£, the magnetic field induces the first-order 
phase transition from B-phase to A- or planar phase at 
the critical field # A B = 0.097rT c0 //i n «3.6kG [U 
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FIG. 2: (Color online) Energy landscape on the unit sphere of 
h, 5Q(h), at (jL n H/nT c0 = 9.2 x 1CT 4 (a) and 0.0061 (b) where 
we fix if/n = —0.5537 which minimizes the dipole interaction. 
We also set H \\ x and T/T c0 = 0.2. The bright (dark) color 
depicts the higher (lower) energy. The energy gap min |£(feii)| 
of Eq. (|3| is displayed in (c). (d) Field dependence of l z 
estimated with the stable (h, if) for Ad /A 2 = 2 x 10~ 4 (the 
solid line) and 2x 10~ 5 (the dashed line). 



Numerical results. — Since the n- vector always points 
to z in the case of the perpendicular field H \\ z. we 
here consider a parallel field H \\ x. Figures and 
[2Kb) describe the energy landscape on the unit sphere 
of h. Figure OJc) depicts the energy gap of the SABS 
evaluated from Eq. ©. The SABS becomes gapless along 
the certain trajectory on the sphere of n, which coincides 
with the condition of l z = 0. The stable configuration of n 
is determined as a consequence of the interplay between 
dipole interaction and Zeeman energy. The former favors 
the situation where n is normal to the surface, namely, 
£ z = 0, while the condition for minimizing the Zeeman 
energy [231 ]. i z = l, is same as the condition that opens 
the maximum energy gap in the SABS. Hence, in the 
magnetic field lower than the dipolar field Hp, m 30G ~ 
0.00l7rT c o/^t n , as displayed in Fig.[5Ja), n points to z . It 
is seen in Fig. [2jb) that for the larger H 's it tends to tilt 
from z to the direction with l z ^0. 

The field dependence of t z estimated with the stable 
configuration of (n,ip) is displayed in Fig. H^d). In the 
limit of the low field, i z is locked to be l z = 0, which 
ensures the existence of surface Majorana fermions. l z 
stays zero up to the critical value [i n H* /itT c q « 0.001, 
which is consistent with the argument that the systems 
with t z = has the discrete symmetry. At H > H*, 
the symmetry protected topological phase with i z = 
undergoes a change to non-topological phase with i z ^0. 

In Fig. EJa) we plot the field dependence of the local 
spin susceptibility on the surface, Xfi Z {z = Q), defined as 
Xttv(z)/XN = M^(z)/M^ for a magnetic field H || f„. 
The local magnetization M M (z) is estimated as M fJi (r) = 
M-^[h t _ l + —^jj (fJn(k, r; iujm))^ u } with that in the normal 



Xzz(0) 



x»(io&)- 



,(0) • 
o 



l; 
l.o 

Z 0.8 
■*i 0.6 
0.4 



0.2 



(a) 



■ m > 

li a H/irT c0 



AL 



(b) 



FIG. 3: (Color online) (a) Field dependence of Xv{ z ) /Xn at 
T = 0.2T c0 . The solid (dashed) lines denote x««(0) (Xzz(W0) 
for H || z and the symbols correspond to Xm^( u ) f° r -H" || a;, 
(b) Af M (z) for if || sb at finH/nT c0 = 9.2 x 10~ 4 (dashed line) 
and 0.0018 (solid lines), where M y , z at ^ n #/7rT c0 = 9.2xlO -4 
are zero. All data are taken with Ad/A 2 = 2 x 10 -4 . 



state Mis 



i , „^ NpH , where Np is the density of states 

ltf 

of the normal 3 He. It is seen from Fig. [3ta) with the 
solid line that for H \\ z 1 the local spin susceptibility on 
the surface, Xzz(0), is considerably enhanced, compared 
with Xz Z (z = 10£) (the dashed line) \3l\ . 

In contrast, when the parallel field (H \\ x) is applied, 
the magnetization M^z) on the surface is sensitive to the 
orientation of £. It is seen in Fig. |3jb) with the dashed 
line that M x (z) at n n H/irT C Q = 9.2 x 10~ 4 is strongly 
suppressed in the surface region, where n \\ £, that is 
l z = 0, is energetically favored. This implies that the 
SABS does not contribute to M x (z) on the surface and is 
consistent with the property of the Majorana Ising spins. 

In the relatively high field ^ n H/nT c ^ = 0.0018, how- 
ever, M x [z) is enhanced around the surface, while M z [z) 
which is perpendicular to H \\ x emerges on the sur- 
face. This emergence of M z (z) on the surface reflects the 
stable configuration of (n,ip), where i z = R xz (n,ip) de- 
viates from zero and is less than unity. As displayed in 
Figs.[3Ua) and[3jb), the magnetic field within the range 
of < £ z < 1 significantly induces M z {z) and x zx {z) on 
the surface, where the SABS opens the finite energy gap 
and the winding number w is not defined. 

Conclusions. — Here, we have clarified that the inter- 
play of dipole interaction and a magnetic field in 3 He-B 
involves a new class of quantum phase transition, that 
is, the topological phase transition with the spontaneous 
breaking of the hidden Z2 symmetry Using the quasi- 
classical theory, we have demonstrated that 3 He-B stays 
topological as the symmetry protected topological order 
and Majorana Ising spins exist unless the Z2 symmetry 
is spontaneously broken at the critical field H* . The 
quantum phase transition is accompanied by the anoma- 
lous behavior of spin susceptibilities, which is observable 
through NMR experiments in a slab geometry. 
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